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1 Introduction 

Big-Bang Nucleosynthesis (BBN) is a pillar of evolutionary cosmology. It is however known 
that "precision" cosmology is not in complete agreement with BBN predictions on fl^h 2 ($7^ : 
baryon density parameter, h : Hubble parameter in units of 100 (km/s)/Mpc ). Agreement 
is however recovered if 7 Li abundance data are disregarded. 

7 Li is the heaviest and most scanty nuclide considered in BBN, but the reaction network 
leading to is hardly questionable, while its abundance data, resumed in the next Section, are 
sound. It seems therefore significant to debate what changes in cosmological parameter 
estimates if 7 Li abundance data are taken on the same foot as other nuclides. 

BBN aims at predicting the abundances of 2 H, 3 He, 4 He, and 7 Li in terms of a single 
parameter: 

m = 10 10 ?i fe /n 7 = 2.7349 x 10 2 n b h 2 . (1.1) 

Here rib an d n 7 are baryon and CMB photons number densities, respectively. Owing to 
baryon conservation and assuming no substantial entropy input between BBN and today, 
their ratio has not changed since then. Data on nuclide abundances will be taken here from 
recent review papers by Matteucci [1] and Iocco et al. [2], while the network of reactions 
used for predictions is accurately illustrated in [2]. 

Tests will be based on WMAP7 data [3]. We shall re-analyse them by using the Cos- 
moMC routines [4] to gauge model likelihood. Together with CMB data, Baryon Acoustic 
Oscillation (BAO) data [5], as well as either Hq [6] or SNIa data [7], as suitably specified, 
will be considered. 

Our analysis will focus on the Dark Energy (DE) state equation, however assumed to 
take a scale independent value w. Neutrino masses will be considered through the parameter 
M u = m u and neglecting their mass hierarchy. A Friedmann-Lemaitre-Robertson- Walker 
metric is assumed, with flat space section. 

2 A quick review on 7 Li stellar abundance 

The history of primordial 7 Li abundance determination starts from the discovery of the 
Spite & Spite plateau [8]. They found that the 7 Li abundance, when the warmest metal- 
poor dwarfs were considered, exhibited no appreciable metallicity dependence, suggesting 



- 1 - 



o 



• Charbonnel & Primas 2005; 'cleanest sample' 
A Asplund et al. 2006 
a Lambert & Reddy 2004 
v Bonifacio et. al. 2007 
O Sbordone et al. 2009 




-1 



HE 1327-2326 
(Aoki et al. 2006) 



-4 



-2 
[Fe/H] 



Figure 1. Observational 7 Li values given by various authors (as in the frame). The long dashed 
black horizontal line is the average 7 Li estimate, taken by [2] and used though this paper. The short 
dashed black horizontal line shows the abundance allowing a best fit among 7 Li and other nuclide 
primeval abundances (see text). The shortest dashed green horizontal lines are the edges of the Ter 
abundance interval consistent with WMAP7 (and related) data. The red and black solid lines show 
the results of galactic evolution models assuming a primeval 7 Li abundance at 2-3<r's below the best 
fit of WMAP7 (and related) data with ACDM models. 

then a primordial value loge( 7 Li) = [ 7 Li/H] = 2.05 ± 0.15. Here [X/H] = 12 + ]og 10 (X/H), 
(X/H) being the ratio of number densities. 

Recent data, using Fe abundances as metallicity indicator, confirm the plateau for 
—3.0 < [Fe/H] < —1.0 but extend to even lower metallicities where a mild further 7 Li 
decrease is detected. This is shown also in Figure 1, extending the outputs from [1]. Ob- 
servational abundances of 7 Li, obtained by [9-13], plus the (controversial) limits in [14] are 
plotted vs. [Fe/H] (l-cr error bars). In the same Figure the results of evolutionary models of 
[15] are shown. They assume a primeval Li abundance as is predicted for a value of t]iq lying 
at ~ 2 — 3 cr's below its best fit for ACDM models, when using WMAP7 and related data. In 
top of it they add the 7 Li contributions deriving from low mass stars, Galactic Cosmic Rays 
(GCR), novae and Asymptotic Giant Branch (AGB) stars (ranging around ~ 41%, ~ 25%, 
~ 9% and ~ 0.5%, respectively). In this way the high metallicity rise observed by [9] is met, 
so confirming that galactic evolution allows a fair understanding of the stellar production 
of 7 Li, but hardly meets its values in the Spite & Spite plateau, let alone lower metallicity 
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Figure 2. Observational abundances of the lightest 4 nuclides produced in BBN, vs. BBN predictions. 
The horizontal cyan dashed bands yield 1-er observational limits. 1-cr WMAP7 (plus BAO and Hq) 
data are comprised between the two vertical dashed lines. The solid vertical line is the best fit with 
all light nuclide data. It meets 3 He and He bands well within 1-er, while it keeps at ~ 1.3 a's from 
both 2 H and 'Li data. 



values. In the plot we added the Li abundance band consistent with WMAP7 and related 
datasets (its 1-a edges are the shortest green lines) and a dashed line representing the average 
of the eight 7 Li estimates quoted in [2] 

( 7 Li/H) = x !0- 10 ; (2.1) 

here the errors are the half-width of the nearly-Gaussian fitting the above eight estimates; 
statistical error on each single estimate are significantly smaller. 



3 BBN predictions on ?] 10 and Q b 

Letting apart cosmological data analysis, we may look at nuclide data only and seek the 
baryon density best approaching a general agreement (g.a.) among all nuclide abundances. 

Within 1-a this is impossible, but, at ~ 1.3 it's from central values, both for 1 Li and 
2 H, a g.a. is met (see Figure 2). This however requires r/io = 5.369, i.e. f]j,/i 2 = 1.963 x 1CT 2 , 
a value significantly below the best fit to WMAP7 and related data. 

Moreover, by looking at Figure 1, we have a visual feeling of the displacement of such 
1.3-cr g.a. point (short dashed black line) from some 7 Li abundance datasets. This point more 
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Figure 3. 65% and 95% likelihood limits on the flbh 2 parameter (derived from WMAP7, BAO and 
Hq data and obtained by fixing w < — 1 values) are compared with the BBN g.a. line (vertical blue 
dashed line) and band (limited by vertical blue solid lines). Decreasing w improves the comparison and 
models with w < —1.8 exhibit an overlap between the 1-cr likelihood interval and the BBN g.a. band, 
including 7 Li. 

or less coincides with observational top estimates of Li abundance and one can hardly escape 
the impression that, in order to agree on this point, one must still invoke some unknown — 
although mild — mechanism for 7 Li depletion in old stars. 

The formal significance of 1.3 cr's is however clear, and tells us that the heavy vertical 
line shown in Figure 2 has a likelihood of some percents. Hereafter we shall refer to the 
short dashed line as "BBN g.a. line". We can also easily determine the interval where 
inside we have the 99% of probability to match all nuclide estimates: 5.25 < rjio < 5.80, 
i.e. 1.92 x 1CT 2 < n b h 2 < 2.12 x 10" 2 ; hereafter we shall refer to this interval as "BBN 
g.a. band". 

The asymmetry of the band is due to the different width of the 1-cr intervals for 2 H and 
7 Li and to the faster rise of 2 H prediction curve in respect to 7 Li. Because of such asymmetry, 
(n b h 2 ) = 2.005 x 10~ 2 . In Figure 2 we also show the 1-cr band (6.031 < rj w < 6.326, i.e. 
2.207 < 10 2 n b h 2 < 2.313), obtainable by applying CosmoMC to WMAP7 data together 
with BAO and Hq constraints, assuming a ACDM cosmology; the plot visually confirms that 
the "BBN g.a. line" lies at ~ 6 cr's from WMAP and related data best-fit value. 

The same plot can be read in a complementary way, by observing that the theoretical 
predictions on 7 Li abundance meet the r/io value best-fitting WMAP7 and related data at 
~ 2.8 cr's (with reference to the data mean and standard deviation shown in eq. 2.1). 

This is consistent with the absence of overlap between the "BBN g.a. band" and the 
ACDM interval; accordingly, the likelihood for 7 Li data and ACDM agreement is <C 0.1%. 
In turn, this shows that agreement would be re-approached if cosmological predictions are 
displaced downwards by ~ 4%. 

In the era of precision cosmology this is a non-negligible shift. However, here we show 
that, if we abandon the assumption that the DE state parameter w = — 1 , a downwards shift 
of such order is obtainable. 
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Figure 4. Distance moduli of Union2 SNIa data are plotted vs. different model predictions. 

In the very WMAP Cosmological Parameters matrix 1 one can easily check that lower 
limits on O^/i 2 can be relaxed, if one extra degree of freedom is opened and w < — 1 is 
allowed. We then performed a number of tests, to determine the best-fit f2t,/i 2 interval with 
fixed w < —1 values, gradually moving away from ACDM. Fits used CosmoMC, constantly 
referring to WMAP7, BAO, and Hq data. Their results are shown in Figure 3. 

Fits were also performed with different assumptions on neutrino mass and we shall 
further comment on this point in the next section. In Figure 3 neutrinos account for a 
fraction f v = 0.07 of dark matter mass. 

Figure 3 shows that decreasing w improves the match between BBN and WMAP7. 
Models with w < —1.8 exhibit an overlap between the BBN g.a. band and the 1-a likelihood 
interval obtained by considering WMAP7, BAO's and Hq data. 

Let us also outline that results are marginally different if SNIa data replace Hq data; 
SNIa and Hq data, however, are not independent and should not be simultaneously used. 
In Figure 4 we however show the distance modulus fx = m — M, obtained for ACDM and 
other cosmologies, against observational distance moduli provided by the Union2 Supernova 
Cosmology Project [7]. The upper solid (magenta) curve, in the main frame, is obtained for 
w = —1.8. Afx's are also plotted, in respect to a ACDM model. There can be scarce doubts 
that a model with w = —1.8 is in perfect agreement with SNIa data. 

4 Phantom w and neutrinos 

There is a wide literature on models with w -C — 1, following the original proposal by [16] of 
an anomalous kinetic energy for the DE scalar field. The main physical context where such 

1 http: / /lambda.gsfc. nasa.gov / product /map / dr4 /parameters, cfm 
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Figure 5. 1- and 2-er likelihood contours on the M v -w plane, when small Vti,h 2 values are set. The 
dotted curves update [18] results, by replacing WMAP5 by WMAP7 data. 

phantom DE was advocated where the limits on neutrino {v) masses (see, e.g., [17]). 

Let then M u = m v be the sum of the i/-mass eigenvalues, for i/'s belonging to the 
usual three particle families. The point we wish to make here is that, if we force f^/i 2 
to shift to values smaller than those obtainable by freely best fitting cosmological data, so 
to approach ' 7 Li constraints, we simultaneously obtain two results: (i) decreasing DE state 
parameters w < — 1 are met; (ii) upper limits on M v are gradually softened. This is shown 
in Figure 5. 

The (i) point is consistent with the results of the previous Section. The only difference 
being that the dataset considered here includes SNIa, instead of Hq data. 

As far as the (ii) point is considered, we show the effects of taking two specific Q^h 2 
values. The former one is the value corresponding to the fit at w = —1.8, discussed in the 
previous Section. The slight discrepancy from the w range found here arises from two reasons: 
(i) Hq data are replaced by SNIa data, (ii) Different parameters are set free: here Qf,h 2 is 
fixed and f u (the fraction of Dark Matter due to z/'s) is free, the opposite choice with respect 
to Figure 3. 

The meaning of the Figure is however clear: If 7 Li data are approached, i/-mass con- 
straint soften. If we opt for the closer model allowing some agreement with 'Li, the upper 
limit on M v shifts from 1.12 eV to 1.27 eV. If we go as down as is needed to meet the g.a. line, 
we find a limit at 1.56 eV. Furthermore, in the latter case, there is an apparent 1.7 cr's "detec- 
tion" for M v . Such pseudo-detection would take a completely different significance if particle 
data would yield a neutrino mass detection in a similar mass range. 
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Figure 6. Rcdshift dependence of the (effective) state parameter in a dynamical DE model. 



5 Discussion 

It may be worth outlining soon that phantom DE is not a straightforward result of quantum 
field theory. If DE is a scalar field, its energy density pde and pressure pde arise from the 
kinetic and potential energy densities (Ek and V, respectively) according to the relations 



p DE = E k + V 



PDE 



E,-V 



(5.1) 



1 



so that w = Pde/ Pde approaches —1 when Ef. <C V, but is unable to bypass the w - 
limit unless E^ lays in the interval —V<E^ < 0. 

Another option is that data analysis tends to favor w < — 1 because Cold Dark Matter 
(CDM) and DE are coupled and the assumed parameter space does not include a coupling 
parameter. In other terms: measuring w <C — 1 could be the consequence of assuming no 
energy flow from CDM to DE if they are actually coupled. 

In this case, data can be approached by using an effective DE energy density [19] 



PDE + Pcdm{1 + zf[f{4>) ~ 1] with C = df/d(t> 
being the CDM-DE coupling (see, e.g., [20]). The state parameter of DE reads then 

,(e//) 



1 



with 



x 



(pcdm/pde)[1 ~ e 



-C(0 o -0)l 



(5.2) 



(5.3) 



(5.4) 



provided that C = (f3/m p ) x 4y / 7r/3 is constant. It must be however outlined that this 
is however an approximate treatment, yielding increasingly rougher results as the coupling 
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increases. In Figure 6 we give an example of the redshift dependence of u/ e ^ f' , for a RP model 
[21] with cosmological parameters close to WMAP7 best fit (Q. c h 2 = 0.114, n h h 2 = 0.0226, 
h = 0.7, A and /3 shown in the frame). In this straightforward case, we find a variable effective 
state parameter which, averaged on z between z = and z = 2.5, yields (w^^) ~ —1.18 . 
As a matter of fact, however, models yielding significantly more negative (w^^) averages 
are hard to build, unless assuming that /3 itself suitably depends on z . The point is that, 
when increasing /3, there arise apparent instabilities in w^ e ^\ even below z = 2 . On the 
contrary, coupled DE models with similar /3's, independently of their capacity to meet data, 
certainly exhibit no such instabilities. Henceforth, the capacity to mimic CDM-DE coupling 
by means of an effective potential is limited. 

However, if w <C —1 is suggested by 7 Li data, this can be generically interpreted as 
possibly favoring CDM-DE coupling. 

In recent work, the option of CDM-DE coupling has been specifically tested against 
data, by simultaneously allowing for higher M u [22, 23], finding that such Mildly Mixed 
Coupled (MMC) models exhibit a likelihood slightly exceeding ACDM (below the 2-a level). 
We plan to reconsider such option in the presence of priors towards lower O^/i 2 . 

6 Conclusions 

The main task of this note is testing the consequence of taking 7 Li data on the same foot 
of other light nuclide abundances. BBN constraints, when 1 Li is disregarded, yield no real 
parameter limitation in top of those arising from CMB, BAO, SNIa, Hq, etc. : the BBN 
interval set by 4 He, 3 He, and 2 H, sometimes set as a prior, more or less overlaps with 
WMAP interval for ACDM cosmologies. 

On the contrary, if 7 Li is considered, BBN is a real extra constraint to cosmological 
models and DE state equations with w <C — 1 are favored together with Q^h 2 values smaller 
by ~ 4-5 %, in respect to the interval usually considered. 

Within this context, M v limits are softened up to 1.6 eV, but the really new feature is 
a quasi-detection for a non-vanishing M v ~ 0.9 ± 0.5 eV. 

A more general, very tentative, conclusion is that there is a realistic possibility that 
WMAP data, for some still unknown reason, led to slightly overestimate f^; if such estimate 
is reduced by 4-5 %, then, a different cosmological scenario seems to open: a simple ACDM 
cosmology is no longer so close to data, predictions of light element abundances by BBN 
would approach a general self-consistency, an energy flow from CDM to DE could mimic a 
unique nature of the dark components, and M v ~ 0.9 eV could turn into an actual detection. 
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